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This  report  describes  a  quantitative  evaluation  of  VLF  and  LF 
propagation  models  used  in  the  WEDCOM  (Weapon  Effects  on  D-region 
COMmunications)  computer  code.  The  evaluation  includes  magnetic 
field  effects,  mode  conversion,  the  effects  of  low  ground  conductivity, 
and  the  effect  of  the  interface  between  VLF  and  LF  propagation  models. 

The  current  WEDCOM  code  gives  good  agreement  with  more  de¬ 
tailed  propagation  calculations  for  many  nuclear  weapon  conditions  of 
interest.  Poor  agreement  with  detailed  calculations  can  result  when 
magnetic  field  effects  are  important,  particularly  when  the  ground 
conductivity  is  also  low.  Mode  'onversion  may  be  significant  for  some 
environments.  Conclusions  and  recommendations  for  model  improve¬ 
ments  are  given. 
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SECTION  1 

INTRODUCTION 


The  WEDCOM  computer  program  (Reference  1)  provides  a  con¬ 
venient  method  for  providing  estimates  of  VLF  and  LF  propagation 
disturbances  that  are  produced  by  nuclear  detonations.  The  distin¬ 
guishing  feature  of  the  program.  1.:.  that  it  combines  an  atmospheric 
chemistry  model,  a  nucleai  ionization  model,  and  propagation  mod¬ 
els  in  a  single  program. 

The  design  of  the  propagation  models  has  been  governed  by  an 
attempt  tc: 

1.  Make  computations  rapidly  so  that  parametric  variations  in 
chemistry  and  ionization  models  can  be  evaluated  economi¬ 
cally  (in  terms  of  computer  usage) 

2.  Achieve  a  balance  between  propagation  computation  accuracy 
and  accuracy  of  the  physical  models 

3.  Allow  a  reasonable  variation  of  inputs  in  terms  of  antenna 
orientation,  ground  conductivity,  path  orientation,  day  or 
night  variations,  etc. 

The  first  version  of  WEDCOM  provided  calculations  for  propa¬ 
gation  between  two  vertically  polarized,  ground-based  antennas. 

The  ground  conductivity  corresponded  to  sea  water  (5  mhos/rreter) 

-2 

or  highly  conducting  earth  (10  mho/meter).  The  nuclear  distur¬ 
bance  was  produced  by  a  single  burst. 
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The  first  major  modification  of  WEDCOM  resulted  in  WEDCOM 
MB.  The  major  additions  were 

1.  Multiple  bursts  (up  to  50)  could  be  included 

2.  Path  conductivity  could  be  given  as  an  input 

3.  Improvements  were  made  in  the  analytic  fit  to  reflection 
coefficients. 

In  response  to  user  inquiries  and  in  anticipation  of  future  uses, 
the  following  changes  have  recently  been  made  to  WEDCOM: 

1.  VLF  propagation  calculations  can  now  be  performed  for  ele¬ 
vated  vertical  or  horizontal  antennas 

2.  An  improved  solution  for  VLF  mode  parameters  when  the 

.3 

path  conductivity  is  very  low  (less  than  10  mho/meter) 
has  been  included. 

The  validation  of  the  code  in  the  past  has  relied  on  the  demon¬ 
stration  that  the  estimates  agree  with  other  more  detailed  calcula¬ 
tions  when  the  nuclear  disturbance  is  intense.  When  the  nuclear 
disturbance  is  weak,  particularly  for  nighttime  conditions  when 
magnetic  field  effects  are  very  important,  the  results  produced  by 
WEDCOM  may  be  questionable. 

The  purpose  of  this  report  is  to  provide  a  quantitative  estimate 
of  the  accuracy  degradation  caused  by  some  of  the  approximations 
and  omissions  in  WEDCOM.  In  particular,  magnetic  field  effects, 
mode  conversion,  and  effects  of  very  low  ground  conductivity  are 
discussed.  The  transition  from  waveguide  mode  theory  to  approx¬ 
imate  ray  theory,  which  is  made  in  the  program  at  a  frequency  of 
30  kHz,  is  also  evaluated. 
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SECTION  2 

CURRENT  WEDCOM  MODELS 


The  documentation  of  the  latest  WEDCOM  propagation  model  is 
being  published  separately  and  the  details  and  equations  will  not  be 
repeatei  here.  Instead,  a  general  description  of  the  model  is  pro¬ 
vided,  followed  by  a  discussion  of.  approximations  that  are  made  in 
the  models.  These  discussions  are  provided  to  introduce  the  com¬ 
parative  calculations  that  are  given  later. 

VLF  MODEL 

The  VLF  propagation  model  is  used  for  frequencies  between  10 
and  30  kHz.  A  waveguide  model  using  impedance  boundary  conditions 
as  described  by  Wait  (Reference  2)  is  used.  The  formulation  is  in 
terms  of  Airy  functions.  Various  limiting  forms  of  the  Airy  functions 
are  used  to  determine  approximate  model  solutions.  The  approxi¬ 
mate  solutions  are  used  as  initial  guesses  to  obtain  a  more  precise 
solution  using  the  Newton-Raphson  method. 

The  impedance  of  the  ionospheric  boundary  is  determined  from 
reflection  coefficients  computed  assuming  a  nonhomogeneous  but  iso¬ 
tropic  ionosphere.  The  ionospheric  impedance  is  defined  to  be  that 
of  an  equivalent  homogeneous  ionosphere  that  would  produce  the  same 
reflected  field  in  the  region  below  the  ionosphere.  An  approximate 
altitude  where  reflection  maximizes  is  defined,  and  the  ionosphere 
boundary  is  assumed  to  be  at  that  altitude.  Analytic  fits  to  reflection 
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coefficients  computed  using  two  real  angles  of  incidence  are  used  in 
an  iterative  solution  to  the  mode  equation. 

Excitation  factors  and  height  gain  functions  for  vertical  dipoles 
are  also  as  defined  by  Wait  (Reference  2).  Excitation  factors  and 
height  gain  functions  for  elevated  horizontal  antennas  are  computed 
using  the  method  defined  by  Galejs  (Reference  3). 

The  electric  field  strength  is  computed  by 

1.  Obtaining  the  isotropic  mode  constant  (attenuation  rate  and 
phase  velocity)  at  various  points  along  the  propagation  path. 

2.  Correcting  the  attenuation  rates  based  on  precomputed  values 
of  attenuation  rates  that  include  effects  of  anisotropy  for 
quasi-TM  modes.  No  correction  factor  is  currently  applied 
for  quasi-TE  modes. 

3.  Averaging  the  phase  velocity  and  corrected  attenuation  rates 
along  the  path. 

4.  Computing  the  field  strength,  using  the  averaged  attenuation 

rates  and  phase  velocities  and  the  geometric  mean  of  excita- 
* 

tion  factors  at  the  transmitter  and  receiver  locations  (com¬ 
monly  referred  to  as  the  WKB  approximation). 

A  vector  sum  and  an  RMS  sum  of  the  field  produced  by  the  individual 
mou,'  <5  .  mputed. 

The  mos  important  approximations  and  omissions  are  discussed 
below. 

Magnetic  Field  Effects 

Anisotropy  produced  by  the  earth’s  magnetic  field  affects  both  the 
excitation  and  attenuation  rate  of  the  waveguide  modes.  Under  iso¬ 
tropic  conditions,  fields  are  linearly  polarized  and  modes  can  be 
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SECTION  2 


i 


labeled  TM  (vertical  E-field  or  TE  (horizontal  E-field).  When  the  ! 

»  ' 

anisotropic  effects  are  moderate,  t,he  fields  are  no  longer  linearly 
polarized,  but  have  a  dominant  component  and  can, be  labeled  quasi- 
TM  or  qua'si-TE.  For  moderate  anisotropy,.  quasi-TE  modes  are 
excited  by  vertical  antennas,  and  quasi-TM  modes  are  excited  by  hori- 

5  | 

zontal  antennas  transmitting  in  the  broadside  direction;  however,  the 

i  ' 

excitation  is  small  relative  to  the  excitation  of  the  dominant  compo-  i 
nent.  Also,  for  moderate  anisotropy,  attenuation  rates  for  TM  modes 
depend  in  a  predictable  way  on  propagation  azimuth.  1 


When  the  ani’sotropic  effects  are  strong,  modes  can  no.  longer  be 
identified  as  quasi-TM  or  quasi-TE,  and  association  of  modes  with 

i  ' 


those  produced  with  isotropic  cases  is  difficult.  , 

In  WEDCOM,  only  modes  that  are  excited  when  the  ionosphere  is 


isotropic  are  considered.  An  attempt  has  been  made  to  associate 


these  isotropic  modes  with  quasi>-TM  or  quasi -T.E  modes  that  occur 

i 

for  a  moderately  anisotropic  ionosphere,  and  to  provide  an  appro¬ 


priate  correction  factor. •  To  determine  the  correction  factor,  exci-  , 


tation  factors  and  attenuation  rates  were  computed  by  L.  A.  Berry, 


Office  of  Telecommunications /Institute  of  Telecommunication  Sci¬ 


ences,  using  the  WEDCOM  normal  ionqsphere  (results  described  in 


Reference  4).  The  calculations  were  performed  with  and  without 
the  geomagnetic  field.  All  calculations  were  performed  for  a  verti¬ 


cal,  ground-based  dipole.  Two  values  of  ground  conductivity  were 


used  (5  mhos/meter  and  10  mho/meter).  For  the  case  including 
the  magnetic  field,  the  calculations  were  performed  for  various  c}ip 
angles,  propagation  azimuths,  and  frequencies.  Determination  of 

1  I  ‘ 

the  necessary  correction  factors  from  the  data  supplied  by  Berry, 
is  discussed  below.  1  1 


i  i 


i 
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Excitation  Factors 

i 

When  all  dip  angleslwere  considered,  thte  magnitudes  of  the  exci- 

•  '  j 

tation  factors  varied  significantly  with  propagation  azimuth.  How¬ 
ever,  by  restricting  interest  to  high  dip  angles  (>45  degrees),  well 

1  '  •'  ,  '  l  . 

excited  modes  could  he  identified  as  quasi-TM  modes  and.thp  varia¬ 
tion  of  excitation  factor  magnitude  with  propagation  azimuth  was 

1  I 

small  (less  than  a  factor  of  2).  Also,  for  the  .well  excited  modes, 

i  '  • 

the  anisotropic  and  isotropic  excitation  factor  magnitudes  were  in!  ■ 
good  agreement.  Therefore,  no  correction  of  the  isotropic  excita¬ 
tion  factors  is  included  in  IvEDCOM.  The  consequences  of  ignoring 

anisotropic  effects  for  weakly  excited  modes  are  more  important  for 

'  2>  ^ 
ground  conductivities  rriuch  less  than  10  mho/meter  and  are'dis- 

•  1  «  i 

cussed  later. 

I 

Attenuation  Rate  Correction  Factdr  i  1 

»  -  , 

Using  the  data  supplied  by  Berry  and  restricting  the  range  of  in¬ 
terest  to  high  dip  angles,  an  analytic  fit  to  the  ratiq  of  attenuation 
rate  with  magnetic  field  to  attenuation  rate  without  magnetic  field 
wa;s  determine^  as  a1  function  of  propagation  azimuth  for  the  first 

t  * 

three  quasi-TM  modes.  .This  ratio  is  the  correction  factor  when 
the  ionosphere  is  undisturbed.  '  , 

i  !  ' 

When  the  nighttime  ionosphere  is  weakly  disturbed,  but  not  de¬ 
pressed  sufficiently  to  cause  anisotropic  effects  to  be  unimportant, 
the  correction  factor  is  assumed  to  vary  exponentially  with  reflection 

i  1  ■ 

altitude.  It  is  applied  ip  full  under  normal  nighttime  conditions  and 
decreased  (or  increased)  tov/ard  unity  exponentially  with  decreasing 

I  v  { 

reflection  altitude.  The  e -folding  distance  for  the  exponential  varia¬ 
tion  was  chosen  to  be  7  km,  because  t,hatiis  approximately  the  e- 
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assumed  that  magnetic  field  effects  are  inversely  proportional  to  the 
ratio  of  electron-neutral  collision  frequency  and  the  electron  gyro- 
frequency. 

No  correction  factor  is  included  for  the  recently  added  quasi-TE 
modes.  A  brief  analysis  of  the  data  computed  by  Berry  (see  Refer¬ 
ence  4)  for  the  quasi-TE  modes  for  high  dip  angles  indicates  that: 

1.  Attenuation  rates  are  not  as  sensitive  to  magnetic  azimuth  as 
quasi-TM  modes 

2.  The  variation  with  azimuth  is  not  regular  enough  to  allow  a 
simple  correction  factor. 

A  mitigating  argument  to  the  crude  treatment  of  the  magnetic  field  is 
that  under  conditions  of  strong  nuclear  disturbance,  the  ionosphere 
is  depressed  sufficiently  that  magnetic  field  effects  are  indeed  un¬ 
important. 

A  discussion  of  the  validity  of  the  approximate  treatment  of  mag¬ 
netic  field  effects  is  given  in  Section  4. 

.i 

Mode  Conversion  *'  - 

Steps  3  and  4  in  the  field  strength  determination  ignore  mode  con¬ 
version  (scattering  of  energy  from  one  mode  into  another  mode)  due 
to  variation  in  ionospheric  properties  (altitude  of  reflection  or  iono¬ 
spheric  impedance)  in  the  direction  of  propagation.  Pappert  has 
shown  (Reference  5)  that  mode  conversion  at  a  day-night  terminator 
is  important  in  determining  field  structure  if  the  transition  distance 
is  much  less  than  1000  km  (distance  to  change  the  reflection  altitude 
by  about  20  km).  Wait  (Reference  6)  has  shown  that  mode  conversion 
is  important  for  changes  in  ionospheric  reflection  altitude  of  a  few 
kilometers. 
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The  importance  of  mode  conversion  for  nuclear-produced  distur¬ 
bances  depends  on  the  weapon  yield  and  detonation  altitude.  In  gen¬ 
eral,  low-altitude  detonations  produce  disturbances  that  vary  rapidly 
along  the  propagation  path,  and  high-altitude  detonations  produce 
widespread  ionization  and  ionospheric  parameters  that  vary  slc.vly 
with  distance. 

An  argument  with  some  validity  is  that 

1.  Small  disturbed  regions  produced  by  low-altitude  detonations 
produce  small  changes  in  field  strength  when  WKB  approxi¬ 
mations  are  used,  due  to  the  limited  extent  of  the  disturbance 
along  the  path.  The  field  disturbance  would  not  be  important 
in  system  analyses  even  if  mode  conversion  is  included. 

2.  The  WKB  approximation  is  nearly  valid  when  the  disturbance 
is  widespread,  resulting  in  a  predicted  disturbance  that  is 
significant;  thus  predictions  that  are  significant  to  systems 
analyses  are  nearly  correct. 

Exceptions  to  the  validity  of  the  WKB  approximation  can  be  found, 
even  for  widespread  disturbances,  and  one  will  be  dicussed  later. 

LF  MODEL 

The  LF  model  has  been  changed  little  from  the  original  version 
(Reference  1).  The  model  uses  geometric-optics  approximations 
with  correction  factors  for  diffraction  effects  and  ionospheric  con¬ 
vergence. 

The  total  electric  field  is  computed  as  the  sum  of  a  ground  wave 
field  and  the  fields  of  sky  waves  with  various  numbers  of  hops.  Two 
field  sums  are  provided:  a  vector  sum  and  an  RMS  sum.  The  latter 
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is  provided  because  of  uncertainty  in  ray  path  length  and  therefore 
phase  when  the  ionosphere  is  irregular  along  the  propagation  path. 

The  LF  model  is  applicable  for  vertical,  ground-based  antennas, 
and  magnetic  field  effects  are  neglected  entirely. 

The  approximate  ray  theory  model  is  readily  adaptable  to  a  non- 
uniform  ionosphere.  A  relatively  simple  iteration  procedure  is  used 
to  determine  the  ray  path  that  connects  transmitter  and  receiver  via 
a  nonuniform  ionosphere.  One  difficulty  is  that  the  geometric -optics 
model  requires  that  a  particular  reflection  altitude  be  chosen.  When 
the  ionosphere  is  intensely  disturbed,  the  reflection  altitude  may  be 
poorly  defined  (i.  e.  ,  significant  reflections  come  from  a  relatively 
thick  region).  Some  calculations  that  show  where  reflections  maxi¬ 
mize  (Reference  7)  indicate  that  the  approximate  definition  of  reflec¬ 
tion  altitude  in  WEDCOM  is  quite  good. 

Results  obtained  using  the  approximate  ray  theory  will  be  com¬ 
pared  with  results  obtained  from  more  accurate  computational  models 
in  Sections  3  and  4. 


SECTION  3 

WEDCOM  EVALUATION 


The  evaluation  has  been  performed  in  two  parts.  In  the  first 
part,  ionosphere  models  have  been  used  that  result  in  reflection  at 
altitudes  low  enough  to  cause  magnetic  field  effects  to  be  very  small. 
With  these  depressed  ionospheres,  the  new  additions  to  WEDCOM 
and  the  approximations  made  in  the  LF  ray  theory  model  can  be 
tested  without  confusing  other  approximations  with  magnetic  field 
effects.  Also,  a  comparison  is  made  of  the  LF  and  VLF  models  at 
the  transition  frequency. 

In  the  second  part,  the  evaluation  has  been  performed  using  the 
WEDCOM  nighttime  ionosphere  model  and  a  hypothetical  nighttime 
ionosphere  model  used  by  Pappert  {Reference  8)  to  evaluate  the 
combined  effects  of  anisotropy  and  low  conductivity  on  vertical  E- 
field  excitation  by  horizontally  polarized  antennas. 

The  electron  density  profiles  that  have  been  used  are  shown  in 
Figure  1.  Profiles  1  and  2  were  used  by  Pappert  (Reference  8). 

Note  that  profile  1  is  similar  to  the  WEDCOM  normal  daytime  pro¬ 
file,  at  least  down  to  the  reflecting  region  near  60  km.  Profile  2 
will  result  in  a  reflection  altitude  significantly  higher  than  would 
result  from  the  WEDCOM  normal  nighttime  profile,  and  thus 
anisotropic  effects  can  be  expected  to  be  larger.  It  is  antici¬ 
pated  that  the  WEDCOM  D-region  chemistry  model  will  be  modi¬ 
fied  in  the  near  future  and  at  that  time  it  may  be  advantageous  to 
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Figure  1.  Electron  density  profiles  used  in  evaluating  WEDCOM 
propagation  models. 


modify  the  nighttime  ionosphere  so  that  an  increase  in  reflection  alt 
tude  results.  Thus,  comparisons  made  using  profile  2  will  provide 
useful  guidelines  for  defining  improvements  or  changes  that  poten¬ 
tially  will  be  required  in  WEDCOM  because  of  the  modified  iono¬ 
sphere. 


VLF  PROPAGATION— ISOTROPIC  CASE 

The  recent  additions  to  WEDCOM  include 

1.  Propagation  from  elevated  horizontal  antennas 

2.  Propagation  over  very-low-conductivity  ground, 

Pappert  (Reference  8)  has  studied  the  effect  of  ground  conduc¬ 
tivity,  ionospheric  anisotropy,  and  antenna  elevation  on  the  excita¬ 
tion  of  vertical  electric  field  from  a  horizontal  electric  dipole.  He 
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considered  a  hypothetical  daytime  ionosphere  and  a  highly  aniso¬ 
tropic  ionosphere  (see  Figure  1).  In  the  notation  of  Wait  and  Spies 
(Reference  9),  the  daytime  ionosphere  hcd  a  /3  of  0.  3  and  a  refer¬ 
ence  altitude  of  70  km. 


The  conductivity  of  the  ground  (a)  was  varied  between  10  and 
10  mhos/meter.  The  relative  dielectric  constant  of  the  ground  (ej 
was  15  for  all  cases,  and  the  magnetic  field  strength  (B^)  was 
0. 44  gauss. 


Pappert  used  a  frequency  of  19.  8  kHz,  a  dip  angle  of  50  degrees, 
and  a  propagation  azimuth  of  270  degrees.  For  profile  1,  reflection 
occurs  at  an  altitude  near  60  km  and  the  results  show  little  effect  of 
anisotropy. 


Figure  2  shows  attenuation  rates  for  the  first  two  quasi -TM  modes 
and  the  first  two  quasi- TE  modes  as  a  function  of  ground  conductivity. 
The  curves  are  from  Pappert's  published  results.  The  symbols  on 
or  near  the  curves  are  from  WEDCOM.  The  results  are  in  quite  good 
agreement.  The  slight  difference  in  the  values  for  mode  3  may  be 
due  to  the  very  simple  fit  used  for  reflection  coefficients  in  WEDCOM 
or  to  slight  effects  of  anisotropy. 

WEDCOM  results  for  the  second  quasi-TM  mode  (mode  3)  were 

_5 

not  obtained  for  conductivities  less  than  4x10  mho/meter.  The 
mode  solution  procedures  in  WEDCOM  do  not  converge  at  the  very 
high  attenuation  rates  associated  with  high  order  modes  and  low  con¬ 
ductivity.  The  attenuation  rate  is  assumed  to  be  equal  to  its  value  at 
-5 

4x10  mho/meter  for  all  lower  conductivities.  Time-consuming 
iterative  procedures  .currently  must  be  used  to  find  the  low-conduc¬ 
tivity  solution.  Because  the  attenuation  rates  are  so  high  when  the 
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Figure  2.  Comparison  of  attenuation  versus  ground  conductivity 
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conductivity  is  very  low,  the  mode  is  essentially  lost.  These  pro¬ 
cedures  are  not  considered  necessary  in  WEDCOM. 

•  I 

Figure  3  shows  excitation  fkctors  for  a  ground-based  vertical 
dipole  as  a  function  of  conductivity.  The  agreement  is  good  for  the 
.quasi-TM  modes.  The  excitation  factor  for  the  quasi-TE  mode  is 

i 

a  direct  result  of  anisotropy  and  cannot  be  produced  by  the  current 
WEDCOM  model.  Note  that  for  this  nearly  isotropic  case,  the 

quasi-TE  modes  are  weakly  excited.  Their  omission  is  of  little  con- 

:  :  I  -4 

sequence  for  conductivites  greater  than  about  10  .  As  pointed  out 

in  Reference  8,  for  very  low  conductivities  the.  weak  excitation  of 

the  quasi-TE  modes  tends 'to  be  balanced  by  the  higher  attenuation 

I 

rates  of  the  quasi -TM  mode,  and  the  quasi -TM  and  quasi-TE  fields 

.  I 

will  be  comparable  for  long  paths. 

Figure  4  shows  the  excitation  of  the  vertical  E-field  versus  ground 

conductivity  for  an  elevated  horizontal  electric  dipole  transmitting  in 
.  I 

the  direction  of  the  dipole  (end  fire).  Again  the  comparisons  between 
WEDCOM  and  Pappert's  result  are  good  for  the  quasi-TM  modes. 

i  'i 

The  excitation  of  the  quasi-TE  modes  is  dependent  directly  on  ani¬ 
sotropy,  and  these  modes 'are  not  computed  for  end-on  launch  in  the 

WEDCOM  model.  The  discussion  given  above  regarding  the  impor- 

.1  ' 

tance  of  the  quasi-TE  modes  applies  here. 

The  results  produced  by  WEDCOM  and  those  described  in  Refer- 

‘  I 

ence  8  for  a  horizontal  dipole  generally  agree  well  when  anisotropic 
effects  are  weak.  However,  thereiare  situations  (low  conductivity, 

1  !  I 

long  paths)  where  even  very  weak  anisotropic  effects  become  im¬ 
portant. 
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Figure  4.  Comparison  of  the  magnitude  of  the  excitation  factor  Aj_|  versus 
conductivity  for  a  horizontal  dipole  at  4  km.  End-on  launch. 
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LF  PROPAGATION-ISOTROPIC  CASE 

The  results  produced  by  the  LF  propagation  model  in  WEDCOM 
have  ^een  compared  to  results  obtained  using  a  computationally  ac¬ 
curate  wavehop  model  supplied  us  by  Berry  of  the  Office  of  Telecom¬ 
munications/Institute  for  Telecommunication  Sciences.  The  wave- 
hop  program,  called  ANIHOP,  is  described  in  Reference  10.  The 
program  requires  as  input  the  four  components  of  the  reflection  co¬ 
efficient  matrix  that  are  computed  by  program  ANIREF,  taking  into 
account  the  anisotropy  produced  by  the  magnetic  field.  ANIREF  is 
also  described  in  Reference  10. 


Field  calculations  as  a  function  of  distance  were  obtained  with 
ANIHOP  using  reflection  coefficients  that  were  calculated  for  propa¬ 
gation  azimuths  of  -90,  -60,  -30,  0,  30,  60,  and  90  degrees  and 
frequencies  of  40  and  100  kHz.  The  static  magnetic  field  was  taken 
to  be  0.  5  gauss.  The  electric  field  strength  variation  with  azimuth 
was  very  small,  as  expected,  because  the  ionosphere  is  nearly  iso¬ 
tropic  for  the  ambient  daytime  ionosphere.  Similar  calculations 
were  made  with  WEDCOM  except  that  no  azimuthal  dependence  is 
obtained  from  WEDCOM. 

Figures  5  and  6  show  the  results  obtained  using  WEDCOM  and 

ANIHOP.  Figure  5  shows  field  strength  as  a  function  of  distance, 

assuming  a  ground  conductivity  of  5  mhos /meter  (sea  water).  The 

WEDCOM  results  are  the  RMS  field  sum.  The  ANIHOP  results  show 

the  vector  sum  of  the  ground  wave  and  wavehops.  Figure  6  shows 

the  same  results  except  that  the  ground  conductivity  had  been  changed 

.3 

to  10  mho/meter. 

The  agreement  between  the  results  at  40  kHz  is  very  good.  The 
WEDCOM  model  predicts  a  significantly  higher  field  than  ANIHOP 
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Figure  5.  Comparison  of  WEDCOM  and  ANIHOP  results  using  WEDCOM 
normal  daytime  ionosphere. 


for  100  kHz.  This  difference  has  been  traced  to  a  difference  in  re¬ 
flection  coefficients,  rather  than  to  propagation  model  approximations. 
The  reflection  coefficients  computed  by  WEDCOM  and  by  AN1REF  are 
compared  as  a  function  of  angle  of  incidence  in  Figure  7.  It  can  be 
seen  that  results  are  in  agreement  at  40  kHz,  but  WEDCOM  predicts 
a  value  of  about  2  dB  higher  than  ANIREF  for  100  kliz.  This  differ¬ 
ence  in  reflection  coefficients  is  sufficient  to  explain  the  difference 
in  calculated  field  strength. 

The  reason  for  the  difference  in  reflection  coefficients  has  not 
been  determined.  It  may  be  caused  by  the  different  numerical  pro¬ 
cedures  used  in  calculating  the  reflection  coefficients,  or  it  may  be 
a  result  of  the  weak  anisotropic  effects  that  remain  while  using  the 
daytime  ionosphere.  Some  calculations  with  ANIREF  using  a  v<*ry 


18 


giT.Ts^yf  jSf^^VTvsr' -i.*  ^  r*5^\  ?,;  .^f^SX-'Tvf  ^T^r  a? 


72TMP-12 


II  l 

Figure  7.  Comparison  of  WEDCOM  and  ANIREF  reflection  coefficients 
using  WE;DCOM  daytime  Ionosphere.  i 

1  '  1  ’ 

low  value  for  the  static  magnetic  field  strength  indicate  that  the 

former  reason  is^  the  most  likely.  ,  ! 

COMPARISON  OF  VLF  AND  LF  MODELS 

AT  THE  TRANSITION  FREQUENCY  >  •'  ; 

Propagation  calculations  made  in  WEDCOM  use  waveguide  mode 

l 

theory  for  frequencies  pf  30  kHz  or  below  and  use  approximate  ray 

:  i  ,  !  * 

theory  for  frequencies  above  30  kHz.  An  isotropic  ionosphere  is 
assumed  in  both  methods.  .  i 

The  waveguide  mode  theory  model  contains  fewer  approximations 

:  i  1  I 

than  the  ray  theory  model;  thus  it  should  provide  a  better  estimate 
of  fielti  strength  than  the  ray  theory  model  when  the  ionosphere  is 

i  ! 

uniform.  On  the  other  hand,  .the  approximate  tray  theory  model  is 

!  i 

more  readily  adaptable  to  calculations  ;with  a  disturbed  ionosphere. 
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Thus  when  the  approximations  in  the  ray  theory  method  are  valid, 
this  method  is  better  than  mode  theory  for  cases  where  the  ionosphere 
is  highly  disturbed.  This  is  mainly  because  mode  conversion  is  not 
included  in  the  waveguide  mode  model. 

The  comparison  between  the  LF  and  VLF  models  in  this  section 
is  for  a  uniform  ionosphere  model.  Calculations  were  made  along 
a  propagation  azimuth  where  the  magnetic  field  correction  factor  is 
unity  in  the  waveguide  mode  calculations.  Comparisons  between  the 
two  models  for  disturbed  cases  are  discussed  later  under  mode  con¬ 
version. 

Figure  8  shows  the  RMS  field  strength  versus  distance  obtained 

at  30  kHz  (VLF  waveguide  model)  and  at  30.  1  kHz  (LF  ray  theory 

model).  Results  were  obtained  for  high  (5  mhos /meter)  and  low 

.3 

(10  mho/meter)  ground  conductivity.  The  ionosphere  used  was 
the  WEDCOM  ambient  daytime  model.  The  agreement  is  very  good 


Figure  8.  Comparison  of  WEDCOM  LF  and  VLF  propagation 
models  using  ambient  daytime  ionosphere. 
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for  high  conductivity  and  good  for  low  conductivity  out  to  about  6 
megameters.  The  difference  beyond  6  megameters  is  due  either 
to  the  use  of  Fresnel  ground  reflection  coefficients  or  to  the  approx¬ 
imation  of  the  antenna  foreground  factors,  both  of  which  depend 
strongly  on  conductivity  in  the  ray  theory  model. 

Figure  9  shows  the  same  comparison,  except  that  the  WEDCOM. 
ambient  nighttime  model  ionosphere  is  used.  Differences  between 
the  two  results  for  lew  conductivity  are  noted  at  3.  5  to  4  mega¬ 
meters. 


In  general,  for  both  ionospheric  conditions,  the  agreement  is  ex¬ 
cellent  considering  the  approximation  made  in  the  ray  theory  model. 


Figure  9.  Comparison  of  WEDCOM  IF  and  VLF  propagation  models 
using  WEDCOM  normal  nighttime  ionosphere. 
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MAGNETIC  FIELD  EFFECTS 


EFFECTS  ON  VLF  PROPAGATION 

The  factor  used  to  correct  mode  attenuation  rates  that  are  calcu¬ 
lated  assuming  isotropic  ionospheric  conditions  was  described  ear¬ 
lier.  The  weaknesses  of  the  method  are 

1.  The  correction  factors  apply  only  for  the  particular  iono¬ 
sphere,  dip  angles,  ground  conductivities,  and  polarization 
for  which  calculations  were  performed 

2.  The  altitude  variation  applied  to  the  correction  factors  has 
not  been  established  by  computations  at  various  altitudes 

3.  It  is  not  always  possible  to  establish  a  one-to-one  corre¬ 
spondence  between  mode  parameters  computed  using  isotropic 
and  anisotropic  ionospheres  (i.  e. ,  the  mode  numbering  sys¬ 
tem  may  be  ambiguous). 

In  analyzing  the  anisotropic  mode  parameter  data  supplied  by 
Berry  (Reference  4)  that  were  computed  using  the  current  VVEDCOM 
nighttime  profile,  some  difficulty  was  experienced  in  choosing  which 
anisotropic  mode  to  compare  with  a  given  isotropic  mode  number. 
However,  in  most  cases,  the  choice  was  obvious,  particularly  at  the 
lower  frequencies  where  quasi-TM  modes  are  well  excited  and  quasi- 
TE  modes  are  poorly  excited  by  a  vertical  antenna.  It  will  be  seen 
later  that  if  the  nighttime  profile  is  raised  a  few  kilometers,  the  ani¬ 
sotropic  effects  become  stronger  and  establishing  mode  correspon¬ 
dence  is  difficult.  Establishment  of  mode  correspondence  becomes 
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even  more  difficult  when  the  ground  conductivity  is  low  or  when  ele¬ 
vated  horizontal  antennas  are  used. 

The  weaknesses  listed  under  items  1  and  2  above  can  be  checked 
and  corrected  by  performing  a  more  complete  set  of  calculations  in¬ 
cluding  anisotropy  that  also  encompass  a  reasonable  variation  of  the 
important  parameters.  These  calculations  would  have  to  be  per¬ 
formed  each  time  the  ambient  nighttime  profile  is  substantially 
changed.  In  performing  the  calculation  to  determine  altitude  varia¬ 
tion  of  correction  factor,  consideration  would  have  to  be  given  to 
the  effects  of  change  in  ionospheric  gradient  as  well  as  altitude.  As 
an  example  of  the  effect  of  varying  ionospheric  altitude,  the  aniso¬ 
tropic  reflection  coefficients  were  computed  for  an  electron  density 
profile  described  by 

N  (h)  -  50  exp  ( 0.  5  (h  -  H  ))  cm  ^  ,  ( 1 ) 

where 

h  =  altitude  (km) 

H  =  reference  altitude  (km) 

R  -3 

Ne<h)  =  electron  density  (cm  )  . 

The  reference  altitude  was  varied  from  60  to  80  km,  which  caused 
the  peak  reflection  altitude  to  vary  over  nearly  the  same  range.  Fig¬ 
ure  10  shows  the  reflection  coefficient  magnitude  of  the  four  compo¬ 
nents  (TEE — vertical  incident,  vertical  reflected;  TEM — vertical 
incident,  horizontal  reflected;  TME — horizontal  incident,  vertical 
reflected;  and  TMM — horizontal  incident,  horizontal  reflected).  A 
frequency  of  20  kHz,  a  dip  angle  of  50  degrees,  a  real  angle  of  inci¬ 
dence  (0.)  of  80  degrees,  and  a  propagation  azimuth  (0  )  of  58  degrees 
1  A 
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Figure  10.  Reflection  coefficient  variation  as  reference  altitude  of 
electron  density  profile  is  varied. 
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were  used.  The  coupling  coefficients,  which  are  a  measure  of  ani¬ 
sotropy,  decrease  with  reference  altitude.  A  curve  decreasing  in 
magnitude  with  decreasing  altitude  with  an  e-folding  distance  of  7  km 
is  also  shown  on  Figure  10.  Based  on  comparison  of  the  TEM  and 
TME  components  and  the  exponential  curve,  the  decrease  of  the  mag¬ 
netic  field  correction  factor  with  altitude  is  too  rapid. 

Figure  11  shows  the  same  plot  but  for  a  propagation  azimuth  of 
238  degrees,  where  coupling  effects  are  stronger.  In  this  case  the 
exponential  curve  is  a  good  representation  of  the  decrease  in  coupling 
coefficients  with  decreasing  reflection  altitude. 

If  the  ambient  nighttime  ionosphere  is  changed  infrequently,  the 
use  of  precomputed  correction  values  may  be  appropriate.  However, 
as  fhe  number  of  variables  is  increased,  resulting  in  increased  dif¬ 
ficulty  of  establishing  mode-to-mode  correspondence,  it  may  be  more 
feasible  to  include  anisotropy  in  the  model  explicitly.  To  minimize 
computer  time,  the  program  could  revert  to  the  isotropic  model  when 
the  ionosphere  is  sufficiently  depressed  to  warrant  the  isotropic  ap¬ 
proximation. 

\ 

As  an  example  of  the  difficulties  in  establishing  mode-to-mode 
correspondence  and  to  estimate  the  effect  of  failing  to  establish  such 
correspondence,  results  obtained  using  the  models  in  WEDCOM  for 
vertical  and  elevated  horizontal  antennas  are  compared  below  to  re¬ 
sults  obtained  by  Pappert  (Reference  8)  for  an  ionospheric  profile 
where  anisotropy  is  very  important.  The  profile  in  the  notation  of 
Wait  and  Spies  (Reference  9)  is  a  f}  =  0.  5,  H*  =  90  profile  (profile 
number  2,  Figure  1). 

In  Reference  8,  a  frequency  of  19.8  kHz,  a  propagation  azimuth 
of  270  degrees,  and  a  dip  angle  of  50  degrees  are  used.  It  should  be 


26 


!  1 


SECTION  4' 


Figure  11.  Reflection  coefficient  variation  as  reference  altitude 
of  electron  density  profile  is  vcried.  • 
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noted  that  the  270-degree  propagation  azimuth  results  in  near  maxi- 

,  i 

mum  coupling  effects  for  the  anisotropic  case. 

The  WEDCOM  isotropic  model  results  were  obtained  using  a  fre¬ 
quency  of  20  kHz.  The  slight  change  in  frequency  should  have  an  in¬ 
consequential  effect  on  the  comparison.  It  should  be  noted  that 
<  >  , 

profile  2  results  in!a  reflection  altitude  (using  the  WEDCOM  defini¬ 
tion)  above  80  km,  which  is  substantially  higher  than  the  altitude  ob¬ 
tained  using  the  current  WEDCOM  ambient  nighttime  profile. 

i  ‘ 

Figure  12  shows  the  attenuation  rates  obtained  (Reference  8)  for 

l 

the  first  five  modes  and  one  high-order  mode.  The  accompanying 
discussion  by  Pappert  indicated  that  there  was  strong  polarization 
mixing  in  all  modes,  and  the  modes  are  not  identified  as  quasi -TM 

*  i 

or  quasi'-TE.  Also  shown  on  the  curve  are  attenuation  rates  for  the 
first  two  TM  modes  obtained  using  She  WEDCOM  model.  Mode  1  has 

J 

the  same  general  variation  with  conductivity  as  Pappert's  modes  2 
and  4,  and  mode  2  has  the  same  general  variation  with  conductivity 

i  ■ 

as  Pappert's  modes  3  and  5.  However,  there  is  no  clear  way  of 

! 

establishing  correspondence  between  the  isotropic  and  anisotropic 
modes. 

! 

Figures  13  and  14  show  the  comparison  between  excitation  factors 
obtained  by  Pappert  and  obtained  using  WEDCOM  models  for  a  ground- 
based  vertical  dipole  and  a  horizontal  dipole  elevated  4  km  and  trans¬ 
mitting  ih  the  end-fire  direction.  Again,  there  is  no  clear  correspon¬ 
dence  between  isotropic  and  anisotropic  values.  Pappert  also 
cotnputed  the  excitation  of  the  vertical  £ -field  from  a  horizontal  an¬ 
tenna  transmitting  in  the  broadside  direction  (not  shown  here).  This 
excitation  factor  does  not  occur  at  all  for  the  isotropic  assumption. 

i 

I  * 

i 

28  i 


ATTENUATION  (d8/Mm) 


SECTION  4 


To  further  examine  the  consequences  of  ignoring  anisotropy  and 
to  see  if  a  method  of  defining  correction  factors  can  be  deduced,  the 
relative  field  strength  as  a  function  of  distance  for  the  various  modes 
was  computed.  Specifically,  was  computed  where 

Edb  =  A  +  8‘  7  ad  '  <2) 

where 
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A  =  excitation  factor  for  vertical  E  (dB) 

a  =  attenuation  rate  (dB/Mm) 

d  =  distance  in  megameters. 

is  related  to  the  actual  electric  field  strength  by  a  constant  and  the 
distance  spreading  factor.  The  constant  and  the  spreading  factor  are  the 
same  for  all  modes  and  for  either  isotropic  or  anisotropic  calculations. 

For  each  type  of  excitation  and  two  values  of  ground  conductivity 
_4 

(1.0  and  10  mho/meter),  E„  was  computed  at  various  distances 

db 

for  all  modes  defined  by  Pappert  and  the  maximum  value  was  deter¬ 
mined.  A  similar  calculation  was  performed  for  the  first  two  TM 
modes  computed  using  VVEDCOM  models. 

Figure  15  shows  a  plot  of  the  values  of  relative  field  strength 
versus  distance  for  a  ground  conductivity  of  1  mho/meter.  The 
curve  labeled  is  the  field  excited  by  a  vertical  antenna  as  com¬ 
puted  in  Reference  8.  The  field  labeled  Vj  is  the  field  excited  by 
a  vertical  antenna  as  computed  by  WEDCOM.  (The  subscripts  A 
and  I  refer  to  anisotropic  and  isotropic,  respectively).  The  num¬ 
bers  by  the  tick  marks  on  the  curves  at  d  =  1  megameter  indicate 
the  dominant  mode  at  short  distances.  The  tick  marks  on  the  curves 
at  other  distances  indicate  where  the  dominant  mode  number  changes, 
and  the  numbers  indicate  the  new  dominant  mode. 

Comparing  the  V  and  V  curves,  it  can  be  seen  that 
A  x 

1.  A  correction  in  both  attenuation  rate  and  excitation  would  be 
required  to  bring  the  fields  into  agreement 

2.  The  change  would  require  that  the  attenuation  rate  computed 
for  isotropic  conditions  be  increased,  contrary  to  the  popular 
notion  that  attenuation  is  increased  for  westerly  propagation 
and  decreased  for  easterly  propagation  as  a  result  of  anisotropy. 
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Figure  15.  Comparison  of  dominant  mode  relative  field  strength  for  the 
vertical  electric  field  at  the  ground  (high  conductivity). 


The  curves  labeled  EF ,  and  EF,  are  relative  vertical  E-fielfis 
i  A  I  !  . 

produced  by  a  horizontal  antenna  at  4-km  altitude  radiating  in  the 

»  ,  1 

end-fire  direction.  A  correction  in  .excitation  factors  'woyld  bring  ( 

these  curves  into  close  agreement.  However,  the  field  excited  from 

the  elevated  horizontal' antenna  in  the  broadside  direction  (labeled 

BS.)  is  greater  than  that  excited  in 'the  end-fire  direction,  and  this1 
A 

field  is  nonexistent  when  the  iono.sphere  is  isotropic. 

i  , 

i 

It  should  be  noted  that  the  mode  association  suggested  by  the  rela¬ 
tive  shape  of  the  attenuation  rate'and  excitation  'factors,  namely  iso¬ 
tropic  mode  2  (2  )  with  anisotropic  mode  5  (5  )  and  1  with  4  ,  . 

,1  A  *  I  A 
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is  in  agreement  with  the1  dominant  modes  for  the  horizontal  end-fire 
antenna,  but  not  for  the  vertical  antenna. 

•  '  f  J 

I  i 

Figure  16  shows  the  same  comparisons  as  described  above  but  for 

■  _4  ■  i 

a  ground  conductivity  of  10  mho/meter,  the  fields  from  the  verti¬ 
cal  antenna  are  in 'better  agreement  out  to  ab<}ut  5  megameters. 
However,  the  change  in  dominant  mode  with  distance  for  the  aniso- 

1  '  1  :  .  i 

tropic  case  would  resiilt  in  field  variations  that  would  not  occur  in 

l  1  1 

the  isotropic  case.  1 

1  ,  ,  , 

The  fields  fromi  the  horizontal  end-fire  antenna  (EF'  and  EF  ) 

A  i  I 

are  in  good  agreement  out  to  about  4,  megameters,;  but  they  change 

*  1  *1 

*  ! 


Figure  16.  Comparison  of  dominant'mode  relative  field  strength  for  the 
vertical  electric  field  at  the  ground  (loW  conductivity). 
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significantly  beyond  that  point.  As  before,  the  vertical  field  from 
the  horizontal  antenna  in  the  broadside  direction  exceeds  the  field 
produced  in  the  end-fire  direction. 

The  general  conclusions  are  that  when  the  anisotropic  effects  are 
as  strong  as  indicated  by  the  examples  chosen  here, 

1.  Any  correction  scheme  would  have  to  include  excitation  fac¬ 
tors  as  well  as  attenuation  rate 

2.  The  choice  of  which  isotropic  mode  to  compare  with  a  given 
anisotropic  mode  to  determine  a  correction  factor  would  be 
extremely  difficult. 

Comparisons  similar  to  the  ones  made  here  will  be  necessary 
when  any  revised  WEDCCM  ambient  nighttime  profile  is  introduced. 
If  the  difficulties  of  defining  correction  factors  for  other  propagation 
azimuths  are  as  great  as  indicated  by  this  study,  it  will  be  necessary 
to  include  explicit  calculations  for  anisotropic  effects  on  WEDCOM. 

EFFECTS  ON  LF  PROPAGATION 

The  electric  field  strength  as  a  function  of  distance  computed  by 
WEDCOM  has  been  compared  to  that  computed  by  AN1HOP  (Refer¬ 
ence  10),  which  includes  the  effects  of  the  geomagnetic  field.  Com¬ 
putations  were  made  with  both  programs  for  frequencies  of  40  and 
100  kHz,  ground  conductivity  of  5  mhos /meter,  and  relative  dielec¬ 
tric  constant  of  80.  The  AN'IHOP  results  were  obtained  for  propaga¬ 
tion  azimuths  of  -90,  -60,  -30,  0,  30,  60,  and  90  degrees.  The 
magnetic  field  strength,  B^,  was  0.5  gauss. 

Figure  17  shows  the  comparison  between  WEDCOM  and  ANIHOP 
computed  fields  for  40  kHz.  The  ANIHCP  results  are  presented  as 
an  envelope  that  contains  the  results  obtained  for  all  propagation 
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azimuths.  The  VVEDCOM  RMS  sum  and  vector  sum  are  shown.  The 
results  compare  favorably  with  the  exception  of  the  2000-to-3000-km 
path  segment.  The  results  are  generally  within  a  factor  of  2  for 
other  path  lengths.  For  this  uniform  ionosphere  case,  the  WEDCOM 
vector  sum  is  to  be  preferred. 

Figure  18  shows  the  field  computed  with  WEDCOM  compared  to 
the  field  computed  with  ANIHOP  for  100  kHz.  Again  the  envelope  of 
the  field  for  azimuths  of  -90,  -60,  -30,  0,  30,  60,  and  90  degrees 
are  given  for  ANIHOP,  and  both  the  RMS  and  vector  field  sums  are 
shown  for  WEDCOM.  For  100  kHz  the  WEDCOM  fields  are  con¬ 
sistently  higher  than  the  field  computed  in  ANIHOP. 

The  reason  for  the  100-kHz  field  difference  is  clearly  shown  in 
Figure  19,  which  shows  the  anisotropic  reflection  coefficient  com¬ 
ponent  magnitude  for  two  propagation  azimuths  as  a  function  of  co¬ 
sine  of  angle  of  incidence.  Also  shown  is  the  isotropic  reflection 
coefficient  computed  in  WEDCOM.  The  comparison  illustrates  that 
the  effect  of  the  magnetic  field  in  addition  to  making  propagation  a 
function  of  azimuth  is  to  generally  reduce  the  reflection  coefficients 
for  all  azimuths. 

Also  shown  in  Figure  19  is  the  reflection  coefficient  obtained  for 
a  very  low  value  of  magnetic  field  strength  in  AN1REF.  (The  value 
used  was  0,01  gauss;  /J^IREF  will  not  run  with  a  zero  magnetic 
field).  This  small  magnetic  field  reflection  coefficient  is  signifi¬ 
cantly  lower  than  the  reflection  coefficient  computed  in  WEDCOM. 

(A  similar  comparison  between  ANIREF  and  WEDCOM  was  noted  in 
Section  3  for  the  ambient  daytime  ionosphere).  It  is  probable  that 
this  difference  is  due  to  different  numerical  procedures  in  the  two 
reflection  coefficient  programs. 
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The  conclusion  to  be  drawn  from  the  results  shown  in  Figures  17, 
18,  and  19  is  that  the  effects  of  the  magnetic  field  can  be  partly  ac¬ 
counted  for,  at  the  higher  LF  frequencies,  by  simply  using  a  reduced 
value  of  the  reflection  coefficient.  The  magnitude  of  the  reduction 
will  depend  on  frequency  and  reflection  height. 


•  I 
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MODE  CONVERSION  ' 
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The  field  strengths  for  each  mode  in  the  waveguide  mode  solution  •  1 

are  obtained  by  using  average  values  of. attenuation  rate  and  phase  1 

velocity  over  the  path  length,  and  using  the  geometric  mean  value  of 

i 

the  excitation  factors  computed  at  the  path  end  points.'  The  appjroxi-  .  > 

mation  is  commonly  referred  to  as  the  WKB  approximation.  1 

.  .  I 

It  has  always  been  clear  that  the  approximation  is  more  nearly  > 
valid  for  some  types  of  nuclear  detonations  (high-altitude,  widespreajd  ! 

debris  sources)  than  for  other  types  (low-altitude,  contained  sources).  , 

When  the  ionosphere  varies  rapidly  in  the  direction  of  propagation,  1  I 

energy  is  scattered  from  one  mode  to  .another  mode.  ’Formulations 
to  compute  mode  conversion  for  certain  idealized  situations  have.  '  ! 

been  developed  (References  6  and  1 1)  and  have  been  used  to  explain 

i  '■  t  I 

the  propagation  across  the  day-night  terminator.  Recently  Pappert  , 

(Reference  5)  has  used  mcde  conversion  to  study  propagation  across  ! 
a  day-night  terminator  as  a  function  of  distance  required  for  the 
transition  from  daytime  conditions  to  nighttime  conditions.  The  re¬ 
flection  altitude  was  assumed  to  vary  from  nighttime  (90  km)  to  day¬ 
time  (70  km)  over  a  distance  AD ,  and,  AD  was  varied  from  0  to  1000 
km.  One  conclusion  was  that  a  transition  distance  near  1000  km  (a 

i 

variation  of  1-km  altitude  per  50-km  path  length)  was  required  before  | 

the  WKB  approximation  and  the  calculations  using  mode  conversion  : 

were  in  substantial  agreement.  1  ! 

i  1 
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Mod$  conversion  is  not  accounted  for  in  WEDCOM,  nor  is  any 
attempt  made  to  recognize  situations  where  the  WKB  approximation 
is  violated. 

In  making  comparisons  between  WEDCOM  LF  and  VLF  models  at 
the  transition  frequency  (30  kHz),  and  for  a  disturbed  ionosphere,  a 

I 

case  was  identified  where  mode  conversion  effects  are  clearly  im¬ 
portant.1  This  case  is  discussed  here. 

! 

The  disturbance  was  produced  by  placing  1  MT  of  fission  debris 
at  an  altitude  of  200  km  directly  over  the  midpoint  of  a  10-Mm  path. 

I 

!The  debris  was  distributed  over  a  315-km  radius,  and  the  computa¬ 
tion  time  was  100  sec  after  detonation.  Nighttime  conditions  were 
assumed  to  prevail  over  the(entire  path. 

I  The  depression  produced  by  this  disturbance  is  shown  plotted  in 
!  Figure  20.  The  reflection  altitude,  as  defined  in  WEDCOM,  is  plot¬ 
ted  as  a  function  of  distance  along  the  path.  Also  shown  is  the  de- 

! 

:pression  profile  4hat  would  b'e  produced  if  the  debris  source  is  moved 
1500;  km  at  right  angles  to  the  propagation  path.  The  steepness  of 
the  sides  of  the  depression  changes  very  little  as  the  source  is  moved, 
although  the  minimum  reflection  altitude  is  significantly  increased. 

I 

Also  shown  for  reference  is  the  nighttime -to-daytime  transition 
analyzed  by  Pappert  (Reference  5)  where  the  transition  distance  is 

1000  km.  The  slope  of  the  transition  region  is  roughly  the  same  as 

i 

that  of  the  nuclear-produced  depression,  but  the  total  altitude  change 
is  much  less.  :  : 

1 

Figure  21  shows  the  field  strength  versus  distance  computed  by 

the  LIf  and  VLF  models  in  WEDCOM  for  frequencies  of  30.  1  and  30 

> 

kHz,  respectively.  Also  shown  for  reference  are  the  fields  computed 
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ionosphere  profile  shown  in  Figure 
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by  the  two  models  for  the  undisturbed  ionosphere.  Despite  the  close 
agreement  for  the  undisturbed  case,  the  agreement  for  the  disturbed 
case  is  poor.  Some  consideration  of  the  mode  str acture  for  the  dis¬ 
turbed  case  leads  to  the  conclusion  that  mode  conversion  may  be  an 
important  contributor  to  the  disparity  in  results.  Those  considera¬ 
tions  are 

1.  In  the  undisturbed  case,  mode  2  is  the  dominant  mode,  even 
to  a  distance  of  10  Mm,  due  to  the  very  weak  excitation  of 
mode  1. 

2.  In  the  disturbed  case,  only  mode  1  can  propagate  any  signifi¬ 
cant  energy  through  the  region  that  is  severely  depressed. 
Thus  mode  2  is  eliminated,  and  much  of  the  difference  be¬ 
tween  disturbed-undisturbed  fields  is  accounted  xor  by  the 
difference  between  mode  2  and  mode  1  excitation  factors. 

If  mode  conversion  causes  significant  energy  to  be  scattered  from 
the  dominant  mode  2  into  mode  1  at  the  beginning  of  the  depression 
(taking  advantage  of  the  lower  attenuation  rate  of  mode  1  in  the  de¬ 
pressed  region)  and  then  if  energy  is  scattered  from  mode  1  back  to 
mode  2  at  the  end  of  the  depression  (taking  advantage  of  the  higher 
excitation  of  mode  2),  the  disturbed  field  strength  may  be  much 
larger  than  the  WKE  value. 

Some  calculations  were  performed  to  see  if  inclusion  of  mode  con¬ 
version  would  improve  the  agreement  between  the  LF  and  VLF  mod¬ 
els.  Mode  parameters  were  available  at  500-km  intervals  along  the 
path  from  the  WEDCOM  calculation.  The  ionosphere  was  modeled 
by  the  series  of  steps  shown  in  Figure  20,  and  mode  conversion  co¬ 
efficients  were  computed  using  the  formulation  described  by  Wait 
(Reference  6).  In  this  formulation  the  field  along  the  vertical  part  of 
the  steps  is  neglected  in  computing  the  conversion  coefficients. 
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The  calculations  started  with  three  modes  incident  at  the  begin¬ 
ning  of  the  depression.  However,  due  to  numerical  difficulties  in 
calculating  height  gain  functions  and  derivatives  of  height  gain  func¬ 
tions  when  the  attenuation  rates  were  large  (resulting  in  large  argu¬ 
ments  for  the  Airy  functions),  some  modes  were  dropped  out  of  the 
calculation  in  the  most  depressed  part  of  the  ionosphere.  The  num¬ 
bers  on  the  steps  show  the  modes  that  were  assumed  to  propagate 
across  the  steps.  Energy  scattered  into  the  higher  order  modes  was 
assumed  to  be  lost  in  the  very  depressed  portion.  Modes  2  and  3 
were  recovered  by  scattering  from  mode  1  at  the  receiving  side  of 
the  depression,  and  mode  2  again  became  the  dominant  mode  at  the 
receiver. 

Figure  21  also  shows  the  field  strength  obtained  as  a  function  of 
distance  for  30  kHz  under  the  disturbed  ionosphere  using  mode  con¬ 
version.  The  agreement  between  the  two  models,  while  still  not 
good,  is  substantially  improved. 

The  case  discussed  here  is  ideal  to  demonstrate  the  importance 
of  mode  conversion,  i.  e. ,  one  where  there  is  a  great  difference  be¬ 
tween  excitation  factors  in  the  undisturbed  part  of  the  path  and  a 
great  difference  between  attenuation  rates  in  the  disturbed  part  of 
the  path.  A  comparison  between  LF  and  VLF  models  for  the  same 
disturbance  but  assuming  the  path  was  in  daylight  showed  much  less 
difference  in  the  models,  because  mode  1  is  well  excited  using  the 
normal  daytime  ionosphere.  However,  this  example  indicates  that 
there  are  situations  where  mode  conversion  can  be  very  important. 

Mode  conversion  can  also  occur  when  the  propagation  path  crosses 
a  land-sea  boundary.  Calculations  performed  by  Wait  (Reference  1-1) 
indicate  that  the  conversion  coefficients  from  one  mode  to  another 
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SECTION  6 

CONCLUSIONS  AND  RECOMMENDATIONS 


CONCLUSIONS  i 

!  tf  I 

i  The  comparison  of  WEDCOM  results  to  results  produced  by  more 
rigorous  models  provides  a  basis  for  some  general  statements  about 
the  computational  accuracy'of  WEDCOM.  The  accuracy  is  strongly 
dependent  on  the  state  of  the  ionosphere  and)  to  a  lesser  extent,  is 

!  >  i  .  : 

dependent  on  ground  conductivity  and  antenna  orientation.  The  con- 
1  .  ,  , 
elusions  are  categorized  below  by  typical  ionospheres: 

1.  ,  Normal  daytime  ionosphere,  or  nighttime  ionosphere  suffi- 
i 

ciently  depressed  to  reduce  anisotropic  effects. 

'  VLF  —  The  WEDCOM  results  generally  agree  well  with  re¬ 
sults  of  more  rigorous  models.  This  statement  applies 

I 

when  the  recently  added  options  that  alloW  calculations -to 

I  ; 

be  performed  for  elevated  horizontal  antennas  and  for 

propagation  over  poorly  conducting  ground  are  exercised. 

An  exception  can  occur  for  long  propagation  paths  over 
*  i  .41 

very  low  conductivity  ground  (<  10  mho/meter). 

i 

LF  —  The  approximate  LF  propagation  results  agree  we. 

I  •  '  ! 

with  the  more  accurate1  mode  calculations  at  the  transition 
■  1  1  1 

frequency  (30  kHz),  and  with)  results  obtained  using  the 

more  accurate  formulation  in  AN1HOP  (Reference  10)  at 

the  .lower  LF  frequency  us.ed  in  the  comparison  (‘40  kHz). 

I  ,  * 

1  j  The  agreement  between  WEDCOM  and  AN1HOP  was  not  as 

,  good  for  100  kHz  as  it  was  at  the  lower  frequencies.  The 
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reason  was  identified  as  a  small  difference  in  reflection 
coefficient  values,  not  to  approximations  in  the  propaga¬ 
tion  model. 

Current  WEDCOM  normal  nighttime  ionosphere 

VLF  —  For  ground-based  vertical  antennas,  the  comparison 
was  made  previously  and  is  essentially  contained  in  Ref¬ 
erence  4,  where  the  derivation  of  the  magnetic  field  cor¬ 
rection  factors  is  presented.  So  long  as  the  program  is 
used  for  vertically  polarized  antennas  and  for  mid-latitude, 
good-to-high  -conductivity  paths,  the  correction  factors  are 
appropriate. 

The  variation  of  the  magnetic  field  correction  factor 
with  decrease  in  reflection  altitude  is  partially  substanti¬ 
ated  by  the  variation  of  the  reflection  coefficient  coupling 
terms  with  altitude  (see  Figures  10  and  11). 

No  results  from  the  rigorous  models  were  available  foi 
comparison  for  horizontal  antennas  and  for  very  low  con¬ 
ductivity  ground. 

LF  —  The  WEDCOM  and  AN1HOP  field  strength  versus  dis¬ 
tance  results  are  in  substantial  agreement  at  40  kHz.  The 
envelope  of  the  AN1HOP  results  (the  minimum  and 
maximum  of  the  field  strength  values  for  various  propaga¬ 
tion  azimuths)  contains  the  WEDCOM  results,  and  the 
maximum  difference  in  field  strength  values  was  generally 
less  than  a  factor  of  2. 

The  WEDCOM  and  ANIHOP  field  strength  results  are 
significantly  different  at  100  kHz.  A  substantial  improve¬ 
ment  in  WEDCOM  results  could  be  obtained  by  simply 
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reducing  the  isotropic  reflection  coefficients  to  account 
for  anisotropic  effects. 

3.  Ionosphere  producing  strong  anisotropic  effects 
(Profile  2  of  Figure  1) 

Only  VLF  was  compared  using  profile  2.  The  reflection 
altitude  obtained  using  this  profile  is  5  to  10  km  higher 
than  that  obtained  using  the  WEDCOM  normal  nighttime 
profile.  Currently  planned  modifications  of  the  WEDCOM 
D-region  chemistry  (Reference  10)  may  be  accompanied 
by  changes  in  the  nighttime  ionosphere  model  that  will  re¬ 
sult  in  reflection  altitudes  between  those  obtained  using 
the  current  nighttime  profile  and  those  obtained  from  pro¬ 
file  2.  Because  of  this  possible  modification,  the  compar¬ 
ison  made  here  indicates  what  modifications  are  necessary 
for  WEDCOM  to  remain  consistent  with  the  new  nighttime 
ionosphere  model. 

The  conclusions  reached  here  are  tentative,  since  the 
comparison  was  made  for  only  one  frequency,  one  dip  angle, 
and  one  propagation  azimuth.  The  tentative  conclusions  are 

a.  Magnetic  field  correction  factors  will  be  difficult  if  not 
impossible  to  define 

b.  Because  of  the  recent  options  allowing  propagation  from 
horizontal  antennas  and  over  poorly  conducting  earth, 
excitation  factors  as  well  as  attenuation  rates  would 
have  to  be  corrected  for  anisotropic  effects. 

4.  Disturbed  profile 

The  comparison  between  results  for  the  LF  and  VLF  models 
at  30  kHz  for  a  severely  disturbed  ionosphere  shows  the 
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potential  significance  of  mode  conversion.  As  described 
in  Section  3,  the  choice  of  disturbance  geometry  and  fre¬ 
quency  tends  to  maximize  conversion  effects.  For  lower 
frequencies  or  for  disturbances  over  path  terminals,  the 
difference  in  excitation  of  modes  is  less,  and  the  effect  of 
mode  conversion  would  be  less  dramatic. 

The  above  comments  apply  to  a  comparison  made  using  two  pro¬ 
pagation  models,  both  of  which  assume  an  isotropic  ionosphere. 
Because  of  the  wide  variation  of  mode  excitation  factors  'j'.d  attenu¬ 
ation  rates  when  anisotropic  effects  are  large,  perhaps  mode  con¬ 
version  will  be  very  important  for  all  frequencies  when  parts  of  the 
path  are  in  the  dark  portion  of  the  earth-ionosphere  waveguide  or 
over  very  low  conductivity  ground. 

RECOMMENDATIONS 

The  comparisons  indicate  that  the  current  WEDCOM  code  gives 
good  agreement  with  more  detailed  propagation  calculations  for  many 
cases  of  interest.  However,  there  are  some  cases  of  practical  con¬ 
cern  where  the  WEDCOM  models  are  deficient  and  improvements 
could  be  made.  The  following  modifications  would  improve  predic¬ 
tions  for  these  cases: 

1.  Add  the  capability  to  compute  reflection  coefficients  taking 
into  account  the  anisotropy  of  the  ionosphere.  This  can  be 
done  by  adapting  the  model  from  ANIREF  (Reference  10)  or 
the  model  used  at  the  Naval  Electronics  Laboratory  Center 
(NELC)  (Reference  12). 

2.  Add  the  capability  to  compute  VLF  mode  constants,  usisg  the 
reflection  coefficients  obtained  by  including  anisotropic  ef¬ 
fects.  Provided  the  modes  can  still  be  identified  as  quasi-TM 
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or  quasi-TE,  the  method  used  in  ANIMOD  (Reference  13)  will 
fit  easily  into  the  framework  of  WEDCOM.  Otherwise,  the 
procedure  used  by  NELC  (Reference  12)  can  be  adopted. 

3.  Retain  the  modified  ray  theory  calculations  for  LF,  but  in¬ 
clude  diffraction  correction  for  elevated  horizontal  antennas. 
The  propagation  model,  given  reflection  coefficients  calcu¬ 
lated  including  anisotropy,  appears  to  be  good  and  is  readily 
adaptable  to  disturbed  ionospheres. 

4.  With  the  more  rigorous  models  included  in  WEDCOM,  retain 
most  of  the  capability  for  fast  computations  by  doing  the  fol¬ 
lowing: 

a.  Organizing  the  program  sc  that  ionization  calculations  are 
made  at  relatively  short  increments  along  the  propagation 
path,  but  so  that  the  points  for  making  propagation  calcu¬ 
lations  are  under  program  control. 

b.  Perform  propagation  parameter  calculations  at  some  ini¬ 
tial  point  on  the  path. 

c.  From  the  ionization  calculations  and  the  path  geometry, 
determine  when  ionosphere  properties  have  changed  suf¬ 
ficiently  to  require  new  propagation  calculations. 

d.  Use  the  mode  solutions  obtained  for  a  previous  calculation 
(the  initial  point  being  the  first  previous  calculation)  as 
inputs  to  obtain  the  solution  iteratively  at  the  next  location. 
The  procedure  will  help  to  establish  mode  correspondence 
in  the  transition  from  anisotropic  to  isotropic  regions. 

e.  When  progressing  from  a  near-normal  nighttime  portion 
of  the  path  to  a  disturbed  region,  revert  to  calculation 
methods  that  assume  an  isotropic  ionosphere.  This  pro¬ 
cedure  should  result  in  performing  a  few  anisotropic 
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calculations  in  the  undisturbed  portion  of  the  waveguide 
and  many  isotropic  calculations  in  the  disturbed  portion 
of  the  waveguide. 

f.  Include  a  rudimentary  form  of  mode  conversion,  or  flag 
results  to  indicate  where  mode  conversion  may  be  im¬ 
portant.  The  art  of  computing  conversion  coefficients 
while  accounting  for  anisotropic  effects,  is  not  developed 
enough  for  WEDCOM  application.  Addition  of  approxi¬ 
mate  mode  conversion  calculations  in  the  isotropic  portion 
of  the  waveguide  probably  would  be  useful. 


53 


BLANK  PAGE 


REFERENCES 


1.  Gambill,  B.  ,  andW.S.  Knapp,  A  FORTRAN  Code  for  the  Calcu¬ 
lation  of  VLF  and  LF  Propagation  in  a  Nuclear  Environment 
(WEDCOM  MB),  DASA  2504,  70TMP-32,  General  Electric—  ' 
TEMPO,  July  1970. 

2.  Wait,  J.  R. ,  Electromagnetic  Waves  in  Stratified  Media,'  Perga- 
mon  Press,  The  MacMillan  Co.,  1952. 

3.  Galejs,  J, ,  "ELF  and  VLF  Fields  of  a  Horizontal  Electric  Di- 
pole,  "  IEEE  Transactions  on  Antennas  and  Propagation,  vol  , 
AP-16,  no.  6,  November  1968. 

:  i 

,  •  i 

4.  Knapp,  W.S.  ,  et  al. ,  Eighth  Annual  Report  on  Development  of 
Aids  for  the  Study  of  High-Altitude  Nuclear  Explosions,  DASA 
2285,  69TMP-18,  General  Electric— TEMPO,  April  1969. 

5.  Pappert,  R.  A. ,  et  al.  ,  Numerical  Results  for  VLF  Mode  Con¬ 
version  in  the  Earth-Ionosphere  Waveguide,  Interim  Report  No. 
712,  Naval  Electronics  Laboratory  Center,  April  1971. 

s 

6.  Wait,  J.  R.  ,  "Mode  Conversion  and  Refraction  Effects  in  the 
Earth-Ionosphere  Waveguide  for  VLF  Radio  Waves,"  J.  Geophys. 
Res.  ,  Space  Physics,  vol  73,  no.  11,  pp  3537-3548,  June  1968. 

7.  Gambill,  B.  ,  and  R.  Rutherford,  WEDCOM  Propagation  Model. 
Improvements,  DASA  2682,  71TMP-27,  General  Electric  — 
TEMPO,  July  1971. 

8.  Pappert,  R.  A.  ,  "Effects  of  Elevation  and  Ground  Conductivity  1 
on  Horizontal  Dipole  Excitation  of  the  Earth-Ionosphere  Wave¬ 
guide, "  Radio  Science,  vol  5,  no.  3,  pp  579-590,  March  1970.. 


i 


55 


72TMP-12 


9.  Wait,  J.  R.  ,  and  K.  P.  Spies,  Characteristics  of  the  Earth- 
Ionosphere  Waveguide  fdr  VLF  Radio  Waves,  NBS  Technical 
Note,  No.-' 300,  November  1.964. 

10.  Berry,  L.A.  ,  and  J.  E.  Herrmann,  A  Wave  Hop  Propagation 

Program  for  an  Anisotropic  Ionosphere,  OT/1TS  RR  11,  U.S. 
Department  of  Commerce,  Office  of  Telecommunications,  In¬ 
stitute  for  Telecommunication  Sciences,  Boulder,  Colorado, 
April  1971.  f 

11.  Wait,  J.  R. ,  "On  Mode  Conversion  of  VLF  Radio  Waves  at  a 
Land-Sea  Boundary,"  IEEE  Transactions  cn  Antennas  and 
Propagation,  vol  AD-17,  no.  2,  pp  216-220,  March  1970. 

I  I 

12.  'Pappert,  R.  A. ,  W.  F.  Moler,  and  L.  R.  Shockly,  A  FORTRAN 

Program  for  Waveguide  Propagation  Which  Allows  for  Both 
Vertical  and  Horizontal  Dipole  Excitation,  Interim  Report  702, 
Navy  Electronics  Laboratory  Center,  San  Diego,  California, 
June  1970.  ' 

I  ' 

13.  Program  ANIMOD — Obtained  from  L.  A.  Berry,  U.S.  Depart¬ 
ment  of  Commerce,  Office  of  Telecommunications,  Institute 
for  Telecommunication  Sciences,  Boulder,  Colorado. 

•  i 

!  ’ 


! 


56  , 


i 


